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. Abstract 

CO 

I A parametrization of the Bose-Einstein correlation function of pairs 

Q ■ of identical pions produced in hadronic e + e~ annihilation is proposed 

within the framework of a model (the r-model) in which space-time 
and momentum space are very strongly correlated. Using information 
from the Bose-Einstein correlations as well as from single-pion spec- 
tra, it is then possible to reconstruct the space-time evolution of pion 
H ! production. 

1 Introduction 

In particle and nuclear physics, intensity interferometry provides a direct 
experimental method for the determination of sizes, shapes and lifetimes of 
particle-emitting sources (for reviews see, e.g., [1-5]). In particular, boson in- 
terferometry provides a powerful tool for the investigation of the space-time 
structure of particle production processes, since Bose-Einstein correlations 
(BEC) of two identical bosons reflect both geometrical and dynamical prop- 
erties of the particle radiating source. Given the point-like nature of the 
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underlying interaction, e + e annihilation provides an ideal environment to 
study these properties in multiparticle production by quark fragmentation. 

2 Bose-Einstein Correlation Function 

The two-particle correlation function of two particles with four-momenta p\ 
and p 2 is given by the ratio of the two-particle number density, p2(pi,P2), to 
the product of the two single-particle number densities, Pi(pi)pi(p2)- Being 
only interested in the correlation R2 due to Bose-Einstein interference, the 
product of single-particle densities is replaced by Po(Pi,P2), the two-particle 
density that would occur in the absence of Bose-Einstein correlations: 



Since the mass of the two identical particles of the pair is fixed to the pion 
mass, the correlation function is defined in six- dimensional momentum space. 
Since Bose-Einstein correlations can be large only at small four-momentum 



difference Q = y—(pi —P2) 2 , they are often parametrized in terms of this 
one-dimensional distance measure. There is no reason, however, to expect 
the hadron source for jet fragmentation to be spherically symmetric. Recent 
investigations, using the Bertsch-Pratt parametrization [6,7], have, in fact, 
found an elongation of the source along the jet axis [8-12] in the longitudinal 
center-of-mass (LCMS) frame [13]. While this effect is well established, the 
elongation is actually only about 20%, which suggests that a parametrization 
in terms of the single variable Q, may be a good approximation. 

There have been indications that the size of the source, as measured us- 



ing BEC, depends on the transverse mass, m t = ym 2 +p 2 = ^ E 2 — p 2 , 
of the pions [12,14,15]. It has been shown [16,17] that such a dependence 
can be understood if the produced pions satisfy, approximately, the (gen- 
eralized) Bjorken- Gottfried condition [18-23], whereby the four-momentum 
of a produced particle and the space-time position at which it is produced 
are linearly related: x = dp. Such a correlation between space-time and 
momentum-energy is also a feature of the Lund string model as incorporated 
in Jetset [24], which is very successful in describing detailed features of the 
hadronic final states of e + e~ annihilation. Recently, experimental support 
for this strong correlation has been found [12]. 

A model which predicts both a Q- and an m t -dependence while incorpo- 
rating the Bjorken-Gottfried condition is the so-called r-model [25]. In this 
article we develop this model further and apply it to the reconstruction of 
the space-time evolution of pion production in e + e~ annihilation. 
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3 BEC in the r model 



In the r-model, it is assumed that the average production point in the over- 
all center-of-mass system, x = (t,r x ,r y ,r z ), of particles with a given four- 
momentum p is given by 

x(p) = arp . (2) 

In the case of two-jet events, a = l/m t where m t is the transverse mass and 

T = \/t 2 ~ * s the longitudinal proper time.* For isotropically distributed 
particle production, the transverse mass is replaced by the mass in the defini- 
tion of a and r is the proper time. In the case of three-jet events the relation 
is more complicated. 

The correlation between coordinate space and momentum space variables 
is described by the distribution of x{p) about its average by 8&(x(p) — x(p)) = 
8&(x — arp). The emission function of the r-model is then given by [25] 

S(x,p) = &tH(t)8&(x - arp)pi(p) , (3) 
Jo 

where H(t) is the (longitudinal) proper-time distribution and pi(p) is the 
experimentally measurable single-particle momentum spectrum, both H{r) 
and p\{p) being normalized to unity. 

The two-pion distribution, p2(pi,P2), is related to S(x,p), in the plane- 
wave approximation, by the Yano-Koonin formula [26]: 

= /dVd'z^pOS^) 1 1 + cos [(p, - p 2 )(x 1 - x 2 )] } . 

(4) 

Assuming that the distribution of x(p) about its average is much narrower 
than the proper-time distribution, Eq. (4) can be evaluated in a saddle-point 
approximation. Approximating the function 5a by a Dirac delta function 
yields the same result. Thus the integral of Eq. (3) becomes 

r&TH(T) P1 (-) , (5) 

Jo Vary 
and the argument of the cosine in Eq. (4) becomes 

(Pi -P2){x\ -x 2 ) = -0.5(airi + a 2 r 2 )Q 2 . (6) 

*Thc terminology 'longitudinal' proper time and 'transverse' mass seems customary 
in the literature even though their definitions arc analogous r = \Jt" — r\ and mt = 
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Substituting Eqs. (5) and (6) in Eq. (4) leads to the following approximation 
of the two-particle Bose-Einstein correlation function: 

R 2 (Q, a u a 2 ) = 1 + ReH fe£\ H M^j , (7) 

where H(u) = J &tH{t) exp^iur) is the Fourier transform of H(t). 
This formula simplifies further if R 2 is measured with the restriction 



cii Pd a 2 



In that case, R 2 becomes 



R 2 (Q,a) = l + ReH 2 i^fj . (9) 

Thus for a given average of a of the two particles, R 2 is found to depend 
only on the invariant relative momentum Q. Further, the model predicts a 
specific dependence on a, which for two-jet events is a specific dependence 
on Wit.' 

Since there is no particle production before the onset of the collision, H{r) 
should be a one-sided distribution. We choose a one-sided Levy distribution, 
which has the characteristic function (Fourier transform) [27] (for a ^ l)* 



H(uj) 



exp 



--(At 

9. V 



OJ 



. an 

- i sign(u;) tan ( — 



+ iujT 



(10) 



where the parameter To is the proper time of the onset of particle production 
and At is a measure of the width of the proper-time distribution. Using this 
characteristic function in Eq. (9) yields 



R 2 (Q,a) 



cos 



ar Q 2 + tan (^pj ^AtQ 



exp 



'aArQ s 



Un the initial formulation of the r- model this dependence was averaged over [25] due 
to the lack of mt dependent data at that time. 
■I-For the special case a = 1, see, e.g., Ref. 28. 
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which for two-jet events is 



R2(Q,m t ) 



1 + cos 



■ exp 



fa-K\ ( ArQ 2 



, tan — 
m t V 2 / \ 2m t 



'ArQ 2 



2m t 



(12) 



We now consider a simplification of the equation obtained by assum- 
ing (a) that particle production starts immediately, i.e., r Q = 0, and (b) 
an average a-dependence, which is implemented in an approximate way by 
defining an effective radius, R = ^/aAr/2, which for 2-jet events becomes 

R = i/Ar/(2m t ). This results in: 



MQ) = 1 ■+ 

where i? a is related to R by 



cos 



(R*Q) 
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exp 



(RQ) 



2a 
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(13) 
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To illustrate that Eq. (13) can provide a reasonable parametrization, we show 
in Fig. 1 a fit of Eq. (13) with R a a free parameter to Z-boson decays generated 
by Pythia [29] with BEC simulated by the BE 32 algorithm [30] as tuned to 
L3 data [31]. In particular, it describes well the dip in i?2 below unity in the 
Q-region 0.5-1.5 GeV, unlike the usual Gaussian or exponential parametriza- 
tions. While generalizations [32] of the Gaussian by an Edgeworth expansion 
and of the exponential by a Laguerre expansion can describe the dip, they 
require more additional parameters than Eq. (13). Recently the L3 Collabo- 
ration has presented preliminary results showing that Eq. (13) describes their 
data on hadronic Z decay [33]. 



4 The emission function of two-jet events 

Within the framework of the r-model, we now show how to reconstruct the 
space-time picture of pion emission. We restrict ourselves to two-jet events 
where we know what a is, namely a = l/m t . The emission function in 
configuration space, S x (x), is the proper time derivative of the integral over 
p of S(x,p), which in the r-model is given by Eq. (3). Approximating by 
a Dirac delta function, we find 

^-™-(?)'*M*(,-=^. < 15 > 
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Figure 1: The Bose-Einstein correlation function R 2 for events generated by 
Pythia. The curve corresponds to a fit of the one-sided Levy parametriza- 
tion, Eq. (13). 

where n and n are the number and average number of pions produced, re- 
spectively. 

Given the symmetry of two-jet events, S x does not depend on the az- 
imuthal angle, and we can write it in cylindrical coordinates as 



where r] is the space-time rapidity. With the strongly correlated phase-space 
of the r-model, rj — y and r = p t r/m t . Consequently, 



where p Puy is the joint single-particle distribution of p t and y. 

The reconstruction of S x is simplified if p Puy can be factorized in the 
product of the single-particle p t and rapidity distributions, i.e., p Pt y = 
Ppt(Pt)Py(y)- Then Eq. (17) becomes 



The transverse part of the emission function is obtained by integrating 
over z as well as azimuthal angle. Pictures of this function evaluated at 
successive times would together form a movie revealing the time evolution of 
particle production in 2-jet events in e + e~ annihilation. 



S x (r,z,t) = P(r, V )H(T) , 



(16) 




(17) 




(18) 
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To summarize: Within the r-model, H(t) is obtained from a fit of Eq. (12) 
to the Bose-Einstein correlation function. From H(j) together with the inclu- 
sive distribution of rapidity and p t , the full emission function in configuration 
space, S x , can then be reconstructed. 
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